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^ Abstract 

^ (— | We consider an exactly soluble model of two Bose-Einstein condensates with a Joseph- 

son-type of coupling. Its equilibrium states are explicitly found showing condensation 
and spontaneously broken gauge symmetry. It is proved that the total number and 
total phase fluctuation operators, as well as the relative number and relative current 
fluctuation operators form both a quantum canonical pair. The exact relation be- 
tween the relative current and phase fluctuation operators is established. Also the 
dynamics of these operators is solved showing the collapse and revival phenomenon. 
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1 Introduction 



Since the 1995-observations |]-[3] of Bose-Einstein condensation (BEC) in trapped alkali 
gases, an intense renewed interest is going on in the research of the physical properties and 
the nature of Bose condensed systems. In particular the interference pattern between two 
overlappping condensates has been measured, see e.g. f| and many other recent experi- 
mental settings and results. 

In this context of interference, the static and dynamic properties of the phase of the 
condensate are of major importance. This has been the subject of many theoretical studies 
all over the last years. As a primordial and old question, the very existence of the phase 
and/or the phase operator, comes into the picture again. 

One encounters continuous efforts to formulate the phase (operator) in the standard theory 
of BEC, which we could call the Bogoliubov-Hartree-Fock theory ||, in a system with a 
finite number of atoms (see e.g. ||). One is constantly assuming that the condensation is 
occuring into a coherent state of the lowest energy mode of the system. Such a state fixes 
a well defined phase and amplitude, but should in stead exhibit inevitable fluctuations 
of both these quantities. Or, one is fixing the number of atoms in the system, i.e. the 
condensation takes place in a number state of Fock space, excluding any atom number 
fluctuations. Although these basic theoretical difficulties are now getting ripe in the minds 
of many researchers in the field, all kinds of procedures and tricks are permanently invented 
to wave away these difficulties. In this paper we take the point of view that these questions 
about the character of the quantum state into which the condensation occurs, and its 
major properties, are nevertheless of major importance. By now it is indeed well known 
that a condensate state is neither a 'pure Fock', nor a 'pure coherent' state in the strict 
mathematical sense, nor in the physical sense. 

As explained above, 'simple' coherent states or Fock states lead to annoying technical 
difficulties in order to describe and understand the essentials of many of the experimental 
challenging measurements on BEC which are constantly performed. After all, condensation 
is up to now, only clearly defined and generally accepted for homogeneous systems. Of 
course, we are aware of different tentatives to introduce decent thermodynamic limits for 
trapped gases. With all this knowledge in mind, we focus our attention here, not on 
the situation of BEC in trapped gases, but on the phenomenon of BEC for homogeneous 
systems, where one has a well defined thermodynamic limit, and where the occurence of 
BEC, accompanied by a spontaneous U(l) - symmetry breaking [0] is well understood. 

Furthermore we take into account that the main entries of the theory of the Bose conden- 
sates and their interference patterns are the particle number fluctuations and the phase 
operator fluctuations. 

The main question is here, can one define rigorously a phase operator fluctuation and 
a particle number operator fluctuation of the condensate? The answer is proved to be 
positive. It is based on the notion of fluctuation operator which was introduced in a 
mathematically rigorous framework some time ago [§,0. We realise however that these 
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results did not reach so far the majority of the theoretical physics community. The aim for 
introducing the notion of fluctuation operator, was precisely to study the quantum effects 
on the level of the fluctuations. We applied this theory of quantum fluctuation operators 
already in order to derive exactly rigorous results on BEC for the Bogoliubov-Hartree-Fock 
model 0. 

In section [|, we describe these results in a language approachable for non-mathematics 
minded readers. The celebrated phase operator is nothing but the canonical fluctuation 
operator conjugate to the number operator fluctuation operator. This section is not just 
for warming up, but it should also shed a different light, than one is used to, on the status 
of the existence and meaning of the phase operator. 

In section [3|, we study a model of two Bose-Einstein condensates with a Josephson-like 
coupling. We look at the static and thermodynamic properties of this model of two con- 
densates. As far as we know, the problem of Josephson oscillations between coupled Bose- 
Einstein condensates has not yet been considered in a mathematically rigorous manner. 
Here we present a solvable model in a full quantum field theoretical setting and in section 
|] and [5] we give a full description and analysis of the dynamical equations of the total 
and relative number operator fluctuations and the total and relative phase operator fluc- 
tuations. About the dynamics we find the exact oscillatory time behaviour of all these 
fluctuation operators. We also detect the so-called collapse and revival phenomena. 

Our work is, as far as we know, the first rigorous one on this topic for homogeneous systems, 
it should also put in a new perspective much of the discussions which are going on in the 



large activity dealing with trapped Bosons and their interference patterns (see e.g. [11-13|). 



2 Number and phase fluctuation operators 

In order to fix our ideas and in view of the model we describe in section ||] for the study of the 
phase interference between two condensates, we present the number and phase fluctuation 
operators for the imperfect Bose gas (or mean field Bose gas) []16,[l7l. We follow the lines 
of [|ITJ but it should be clear that its validity is much larger JTJ]. 

The leading idea of this section is to make clear that the up to now rather 'mysterious' 
phase operator, which everybody uses in the field, but about which there are doubts on its 
very existence, has a firm mathematical definition in an equilibrium condensed state of a 
Bose gas. It should be realised that such a condensed state is neither a coherent, nor a Fock 
state in the technical narrow sense. It is defined as the fluctuation operator canonically 
adjoint to the number fluctuation operator. This definition is not just a formal thing, but 
it gives a physical interpretation of the phase operator, different from the existing ones. 
We are not repeating here all the mathematics of the definition of the phase fluctuation 
operator, which can be found in various papers (see e.g. P,|T5|]). We content ourself here 
in making these definitions plausible for the imperfect Bose gas. 
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Let A C I 3 be the centered cubic box of length L, with periodic boundary conditions. The 
Boson creation and annihilation operators in the one-particle state ipL,k( x ) — V~ l ^ 2 e %k ' x , 
x G A, k G A* = ?fZ u are given by 



-ik.x 



a* Lk = — = I dx a*(x)e tk ' x and az,k = —= f dx a(x)e 
yV J a VV J a 

with 

[a(x),a*(y)] = 5{x - y). 
The imperfect Bose gas is specified by the local Hamiltonian Hl ]Tj|: 

H L = T L - f, L N L + ^Nl (1) 

where 

J-L — e k a L,k a L,k , £k — 2^ 

/cGA* 

N L = a L,k a L,k] 

fceA* 

A > measures the strength of the mean field inter particle repulsion. 

This model is exactly soluble in the thermodynamic limit L —>■ oo, keeping the particle 
density in the Gibbs state u)l(-) for ([]]) constant, equal to p, i.e. for all L: 

^u L (N L ) = p. 



It is proved rigorously [ITS] that for T < T c or p large enough, the limit state up{ 
hrn^oo exists as an integral over ergodic states (u/g, a G [0,27r]): 



with 



where 



1 f 27r 
27r Jo 



(jJ a^[a*{f)+aU)]^ = e ~^(f,K f)+2i^\f (0)\ cos a ^ 



1 2,,.\ ..xl ^ * ^ r2nau\ 



Kf(k) = -f(k)coth^, fEL 
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The spontaneous gauge symmetry breaking accompanying the phase transition is visible 
in the states having the property: 



lim ^(^) = v /p^' 



Clearly, p is the condensate and a is the phase of the order parameter. One also proves 
that in the state uig, one has the operator limit: 



a} n 



From now on we limit our attention to one of the ergodic states Up for some fixed a, and 
without restriction of generality we take a = 0, and with a condensate density p ^ 0. For 
simplicity, denote the state tv® by u>. 

The state cu does not have the gauge symmetry. The generator of the gauge symmetry is 
the number operator 

Nr= I dxa*(x)a(x) 



A 



with local number density operator n(x) = a*(x)a(x). The common choice of order pa- 
rameter operator is V^ 1 ^ 2 J , LQ , or taking a self-adjoint combination: 

L = — =(a* L0 - a L . ) = — / dx (a*(x) - a(x)), 

with local order parameter density operator o{x) = i(a*(x) — a(x)). 

We concentrate now on the k-mode fluctuations, with k ^ 0, of the local number and order 
parameter density operators, i.e. on 

FL,k{ n ) = —/= I dx — uj(n(x))) cosk.x 



FL,k{o) = —i= I dx (o(x) — uj(o(x))) cosk.x. 



A 

Remark that for all finite L, the quantities Fi^(n) and Fl^{o) are operators and do 
represent the fluctuations of the number density and of the order parameter density. 

The first tedious question that is posed, is to characterize the limit operators: 

F k (n) = lim F L)k (n) 
F k {o) = lim F L>k {o). 

L— >oo 



The details of the proof of these limits can be found in [10]. Here we just mention that 



the limits are taken in the sense of a central limit theorem. The main result is that the 
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limits F k {n) and F k (o) are operators on a well specified Hilbert space, H k , generated by 
a normalised vector Q k and vectors F k (Ai) ■ ■ ■ F k (A n )Q k , with the A* local operators, like 
e.g. n(x) and o(x), and with arbitrary n; the scalar product of 7i k is given by 

(Ffc(Ai) • • ■ F fc (A„)fi fc , F fc (Bi) ■ ■ ■ F k (B n )n k ^ = 5„, m Perm [ (n k , F k {Ai)F k (Bj)Vt k ^ ] . 

^(2) 

with two-point function given by 

(n k , F k (A i )F k (B j )Q k ) = lim u){F Ljk {Ai)F L)k {Bj)) (3) 

i.e. essentially determined by the two-point functions of the given state u. Therefore as 
is clear from (0), all (n + m)-point functions are given by the two-point function (|3[). The 
definition (0) defines completely the fluctuation operators F k (A) on the Hilbert space 7i k . 

On the other hand, the non-commutative law of large numbers, here of large operators, 
leads straightforwardly to the canonical commutation relation 

lim \F L k (n),F L k (o)] = lim — / dx [n(x), o(x)) = co([n(x), o(x))) = i^, 

L—ioo L^oo Z V J ^ 

or 

[F k (n),F k (o)]=i^. (4) 

This is the basic result for the definition of the phase operator. Equation (£|) means that 
the number fluctuation operator F k (n) and the order parameter fluctuation operator F k (o) 
are canonically conjugate (compare with [q,p] = ih). Hence for the physics of BEC we 
found on the level of fluctuations the canonical pair (F k (n) , F k (o)) . Clearly the operator 
F k (o) satisfies all basic physical requirements for playing the role of what is usually called 
the phase operator of the condensate. 

The reader will have recognised from (^|) and (]4]) that the fluctuation operators 

F k (A), F k (B), . . . 

form an algebra of Boson field operators and in particular that F k (n) and F k (o) form 
quantum canonical variables with a quantisation parameter ^J~p^ (compare with K). On 
the other hand, from (Q) it is clear that the vector fl k defines a generalised or quasi free 
(gaussian) state 

on the Boson field algebra (see [[H|j ). 

This means that the central limit theorem for the k-mode fluctuations in the state u defines 
an equilibrium state uj k on the fluctuation operators. The mentioned quasi free character 
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means that all correlation functions of limit fluctuation operators are polynomial functions 
only of the one- and two-point functions. For more details we refer once more to HlQfl . 

The reader might ask for the unicity of this phase operator, if being defined only as the 
canonically adjoint operator to the number fluctuation operator. The interested reader is 
referred to section and jl0|,|r4| for that discussion. 



3 The model and equilibrium states 

We consider two coupled Bose-Einstein condensates, each of them modelled by an imperfect 
or mean field Bose gas. Denote af(x), i = 1,2 the creation and annihilation operators for 
the two Bose gases, i.e. 

[ai(x),a*(y)] = 5 itj S(x - y). 

We assume that the two gases have the same particle density | (hence p is the total particle 
density p = ^ujl(Nl) with Nl = N\ t L + N^l), and also that they are of the same type of 
particles (i.e. there is only one mean field constant A). We also assume a phase difference 
(p between the gases, and model the Josephson coupling between the gases by a term 

= -7 a tk^, k e- iip + o^ai,*^, (5) 
fceA* 

with 7 > the coupling constant. 

Hence the local Hamiltonian of the system we study is given by: 
H L = T ljL + T 2iL - p L N L + -^N 2 L + C l L 2 

~ 7 Yl a *i,k a ^e~ lLp + a* 2>k a 1>k e lv . (6) 
fceA* 

In this section we find the limiting Gibbs states u = limi^oo ujl at inverse temperature /3 
of this model. A rigorous study along the lines of |18| is perfectly possible, but we permit 
ourselves here a more intuitive approach. As in any mean field model, we replace the 
Hamiltonian (||) by a state dependent effective Hamiltonian 

H t= Yl ( £fc ~ ^ + X P)( a tk a i,k + a*2,k a 2,k) ~ 7 Y a h a 2,ke~ llp + a\ k a 1M e lv (7) 
fceA* fceA* 

with p = linif^oo p L in correspondence with the constraint p = yu(N L ). This effective 
Hamiltonian is bilinear in the creation and annihilation operators and therefore it can be 
diagonalized by a Bogoliubov transformation. 
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Let = [H%, ■} be the generator of this dynamics and f k = £k — fJ> + Ap. Then 



st [ i' k ) = [ T )[ i' k 



The matrix 



fk 

--ye-* f k 

has eigenvalues E^, 

E± = f k ± 1 , 
with corresponding eigenoperators b*± k , 



^ = i(^ T ^), (8) 

i.e. 

Wi >fe ) = ^63=,*- (9) 

The still satisfy Boson commutation rules. The energy spectrum of the quasi-particles 
b±(x) has two branches, {E^\k el"}. 

By a standard argument, using the Boson commutation rules and the correlation inequal- 
ities [pO|,[n[|, characterizing the limit equilibrium states, 

lim (3u(X*5" L (X)) > u(X*X)\n^^ (10) 
i^oo v ' uj(XX*) 

for X any polynomial in the creation and annihilation operators, one finds in a straight- 
forward manner: 

u(b*± lk b±,k) = 8E l • (11) 
e pE k - 1 

Along the usual lines of the derivation of Bose-Einstein condensation, we find a critical 
density (or inverse temperature) above which one derives the following value of the chemical 
potential: 

P = Ap - 7, 

i.e. E^ = efc and E^ > 2j for all k. Hence there is a macroscopic occupation of the 
0- momentum state of the '—'-mode possible; the condensate density is given by: 

lim -u;(&* „&_,<,) = Po >0. (12) 

V— >oo V 
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There is no condensation for the '+'-mode since lim^o E£ = 27 > for \x = Xp — 7. 
From p = lim^oc V~ 1 uj(Nl) and @, one finds 

p= lim ^J2 u ( b *-k b -k) +^(K,k b +,k) 

L — >oo y — 

k 

f dk 1 f dk 1 

= Po+ 77^7^7 7 + 



For a > 0, denote 



then (|T2"|) becomes 



„ (2vr) 3 - 1 y RB (2tt) 3 e^+^i) - 1 ' 



P(a)= C (2tt)3^)-1' (13) 



Po = p-p(0)-p(2 7 ). 
It is clear that p(0) + ^(27) is the critical density. 

The Bose-Einstein condensation ( [12]) implies a spontaneous breaking of the gauge symme- 
try, i.e. the limiting Gibbs state uj decomposes with respect to the U(l) gauge group into 
distinct extremal equilibrium states: 

1 f 27T 
uj = — / dO u>0, 
2?r Jo 

where each of the states ujq is determined by the two-point function (|TT|) and the one-point 
function 

Lo e (bl(x)) = y/p- e ie . (14) 

This important point, whose non-triviality is often overlooked, is proved in J7|. Of course, 
the other way round, namely that gauge symmetry breaking implies condensation, is well 
known and follows trivially from the Schwartz inequality. 

From now on we choose a particular extremal equilibrium state cog, and without loss of 
generality we take 9 = 0. For notational simplicity, this state is again denoted by uj. Using 
once more the correlation inequalities fli~0|), it is not difficult to show that the higher order 
correlations decompose into sums and products of one- and two-point correlations, given 
by (|l|D respectively (0), i.e. the state uj is quasi-free. Therefore we have completely 
characterized the equilibrium states. 

It is clear that the gauge symmetry breaking state under discussion indeed corresponds to 
a state of two different condensates interacting through a Josephson coupling (^). Since 
uj(b* + (x)) = and uj(b*_(x)) = y/po, we find using @ 



uj(al(x))e—^ = uj(a* 2 (x))e^ = J^. (15) 
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Notice that ip is indeed the phase difference between the condensates. Our arbitrary choice 
9 = then actually determines both phases to be ±^ and the choice of equal particle 
densities | yields equal condensate densities y. Moreover it is obvious that the gap less 
mode is related to the broken gauge symmetry, i.e. to the Bose-Einstein condensation, 
and that the mode E£ with energy gap 27 arises due to the presence of the Josephson 
coupling. The detailed study of the fluctuation operators corresponding to these two 
excitation branches and modes is the subject of the subsequent sections. 



4 Total number and phase fluctuation operators 

Motivated by the discussion in section and ||10|| , define the total number and phase 
fluctuations in the box A by (k 7^ 0) 

FL,k(n tot ) = —j= I dx(a\(x)ai(x) + a* 2 (x)a2{x) — p) cosk.x (16) 



W J A 

F L ,k{4>tot) = -7= / dx(b*_(x) - b-(x)) cosk.x, (17) 

V V J A 

where b-(x) is defined in (H). 

Again on the basis of the law of large numbers: 

lim [F Ltk (n tot ) , F L<k ((j) tot )] = lim — = (&* + 6_ ) = iy/po- (18) 



Using (H) one writes also 

a\{x)a\{x) + a* 2 (x)a2(x) = b*_ + b* + (x)b+(x), 



or 



F L ,k(n tot ) = F L)k (n^) + F L:k (n +) . 



with n± defined in the obvious sense. Hence the total number operator fluctuation is 
the sum of two number operator fluctuations of two imperfect Bose gases. For a single 



imperfect Bose gas (see section |2|), one finds this analysis as the subject of [1C], making 
the present analysis straightforward. 

As already discussed in section |^, and apparent from (|T8"D, the limiting fluctuation operators 
F k (-) satisfy Bosonic commutation rules (although the local fluctuation operators do not). 
Equation ( |TB"D learns also that the fluctuation operators F k {n tot ) and F k ((j) tot ) constitute a 
canonical pair, generating an algebra of canonical commutation relations (CCR) of fluctu- 
ation observables of the system. Furthermore, it is shown in and briefly discussed in 
section 13 that the central limit theorem also fixes an equilibrium state u> k on this algebra 
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of limiting fluctuation operators, which is a CCfl-algebra of Bosonic field operators. This 
state is shown to be quasi-free and gauge invariant, and hence completely determined by 
its two-point function, given by HQ: 

u k (F k (A)F k (B)) = lim / dz u(A(z)B(0)) cos k.z, 

L^ccJ A 

where A, B are (in the present case) polynomials in the microscopic canonical Bosonic 
field operators. Remark that there are no technical problems related to the central limit 
theorem for k 7^ 0, off-diagonal long-range order correlations do appear only at k — 0. 

The first step in our study of the total number- and phase fluctuation operators is to 
determine their variances. 

Proposition 1. For k 7^ we have 

(i) LO k (F k (n tot ) 2 ) = u k (F k (n_) 2 ) + u k (F k (n + ) 2 ) 

p Be k If dp e Pe v+ k + e? e * 

coth — — h 



2 2 2 V (2t^) 3 (e**+* - l) (e^ - l) (19) 

1 r dp e ^K+k + e ^+ 
+ 2 ] v {2nf - 1) ( e /3< - 1) 

(ii) Co k {F k {<p tot f) =iooth^. (20) 

Proof, (i) The fact that 6j k (F k (n to t) 2 ) = ^k(F k (n ) 2 ) + u k (F k (n + ) 2 ) follows from the fact 
that the '+'- and '—'-mode are independent of each other. For a calculation of the explicit 
expression for cu k (F k (n to t) 2 ) , see |10| . 



(ii) This is simply the two-point function of the state u, see [jT0" |. □ 



Remark that the result of this proposition, supplemented with ( Jig ) is sufficient in order to 
characterize completely the limiting fluctuation operators F k (ntot) and F k ((p to t) on a well 
defined Hilbert space (see section |2|). We do not enter into these technical details. 

Of course we are interested particularly in the limit k — > of these operators, in order to 
see how the long-range order due to the Bose-Einstein condensation manifests itself on the 
level of the fluctuations. In [ID] one can find a discussion demonstrating that quantum 



effects will only be present in the limit A; — > if one works in the groundstate u 9 , defined 
as the zero-temperature limit of the equilibrium state u: 



uj 9 = lim u. 



At T = 0, the variances fll9|) and (^0|) simplify to 



Po 
2 

ul(F k (Kt) 2 ) = \, 



u 9 k {F k (n tot ) 2 ) 
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and the limit k — > is trivial: 

F (n tot ) = Urn F k (n to t) and F ((p tot ) = lim F k ((p tot ). 

k^O k^O 

These are well defined fluctuation operators, satisfying 

[Fo(ra tot ),Fo(0 to t)] =iy/po 
ti 9 (F (n tot ) 2 ) = f 

^(F (^) 2 )=^. 

We derived the exact uncertainty relation between the number operator and phase operator, 
given by 

^(^K 4 ) 2 )^(F o (0 tot ) 2 ) = ^. 

Remark that the condensate density po, in fact ^/po, acts in this equation, as well as in 
equation (|T8|), as a quantisation parameter (compare with h). Consequently, the whole 
content of our results, as well as all physical interpretations, do disappear in the absence 
of condensation, i.e. if p = 0. 

Finally remark that we have omitted in our notation the state dependence of the fluctua- 
tion operators throughout, although this dependence is important. Fluctuation operators 
corresponding to different states (e.g. corresponding to different temperature or different 
phase) in fact are not comparable as they act on completely different Hilbert spaces. We 
do not enter into these mathematical subtleties. 



5 Relative number and phase fluctuation operators 

The relative number operator in a finite volume is: 

N L ,rel = N liL - N 2iL , 

and its fc-mode fluctuation 

Fl k{n re i) = 7777T f dx (a*Ax)ai{x) — an(x)a 2 (x)) cosk.x. 

V 1 J A 

As before, we are primarily interested in the limit L — > oo, followed by the limit k — > 
0. The relative number operator Ni^ei is not the generator of a symmetry of the local 
Hamiltonian Hl (§) because of the Josephson coupling term. This means that there is no 
question of spontaneous symmetry breaking for the relative number operator, as it is not 
a symmetry. Furthermore a straightforward computation of the dynamics of Ni^ei learns 
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that its spectrum belongs to the excitation branch E + — E~ . Since E£ — = 27 > 0, 
the spectrum of N L ^ re i shows an energy gap. Hence the 0-mode fluctuations of the relative 
number operator and its adjoint will be normal. Therefore the limits L —>■ 00 and k —>■ 
may be interchanged and hence the starting point of our investigation can be the operator 
(i.e. the case k = 0): 

F L {n r ei) = / dx (a*(x)ai(x) - a* 2 (x)a 2 (x)) = — ^ £ a\,k a i,k ~ a\k a 2,k- (21) 

•^ A k 

Next define the fluctuation operator of the relative current: 



F L (jrei) = 27 yi/2 / dx {a\(x)a 2 {x)e lip - a%{x)ax{x)e ltp ) 



A 



T12 a Xk a 2,k e t<P ~ a*2,k a i,ke lip ■ (22) 



27 V 1 / 

This operator clearly corresponds to the 0-mode fluctuations of the relative current from 
one gas into the other. Below, we show in a series of steps that the operators (|21|) and ( p2|) 
are each others adjoint in the limit L — > 00. Afterwards we show the relations between the 
relative current fluctuation operator on the one hand and the relative phase fluctuation 
operator on the other hand. 

A central limit theorem and reconstruction theorem can again be proved for these operators 
(see ||[1(J), proving the existence of the Bosonic field operators 

F{n re i) = lim F L (n re i) 

L^oo 

F(jrei) = lim F L (j rel ) 

L^oo 

in a rigorous mathematical sense. 

Let <L(-) = hniL^oo [H%, •] be the infinitesimal generator of the dynamics (H^ is the effective 
Hamiltonian defined in equation (0), with fi = Xp — 7). 

In order to prove the properties below, it is convenient to write ( pT|) and (^) in terms of 
the quasi-particle operators (§): 

F L (n rel ) = 77^ £ 6 + + h -,kKk (23) 



k 



F L (jrei) = £ b* +tk b^ k - b*_ >k b +ik . (24) 



k 



Proposition 2. The operators F(n re i) and F(j re i) form a canonical pair and satisfy 

[F(n re i),F(j re i)] = ic re i, 
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where 



c re i = (3 lim (F L (n rel ), F L (n rel )) = -(p + p(0) - ^(27)) > 0; 

L— >oo v 7 7 v ' 



p(a) is defined in ^Wj), and (-, zs i/ie Duhamel two-point function defined below. 



Proof. The first statement follows again from the general theory on normal fluctuation 
operators |9, Kj. Also, it is an easy calculation to show that 

[H£,F L (j rel )] =iF L (n rel ). 
Because of the presence of the energy gap 27 > 0, this can be written as: 

F L (jrei) = iS~ l (F L (n rd )) . 
Therefore, using the fact that oj is an equilibrium (KMS) state, one gets 
\F(n re i),F(j rel )] = lim u( \F L (n re i), F L (j r ei)] ) 

= 1 lim u(F L {n rel )[l - e-^]5j 1 F L (r2 rei )) 

= i(3 lim (F L (n re i), F L (n rel )) . 

The explicit expression for the Duhamel two-point function (i^(n re /), FL(n re i)) then fol- 
lows from 

u([F L (n rel ),F L (j rel )]) = JLj2u(b*_ >k b„, k - b* +>k b + , k ) ^> -( Po + p(0) - p(2 7 )). 

^ k ^ 

□ 

The infinitesimal generator 5 U of the microdynamics induces a natural infinitesimal gener- 
ator 5 W of a dynamics on the macroscopic fluctuation operators by the formula : 

5 U F(A) = F(S U (A)). 

Proposition 3. The infinitesimal generator 5 W on the macroscopic fluctuations is given 
by: 

LF(n rel ) = -z(27) 2 F( Jrei ) (25) 
~^F{j rel )=iF{n Tel ). (26) 

Hence F(n re i) and F(j re i) are eigenvectors of <5 2 : 

~8lF{n rel ) = (2 7 ) 2 F(n re/ ) , 5%F(j rel ) = (2 7 ) 2 F( Jre/ ), 



14 



yielding the macrodynamics a>t on the fluctuation operators: 

a t F(n rel ) = e itL F(n rel ) = F{n rel ) cos(2 7 t) + (2j)F(j rel ) sin(2 7 t) 

& t F(j rel ) = e itL F( 3rel ) = ^-F(n rel ) sin(2 7 t) + F(j rel ) cos(2 7 t). 

2 7 

Proof. This follows immediately from the relations 

[HZ,F L (J rel )] =iF L {n rel ) 

[H%, F L (n re i)] = -^(2 7 ) 2 F L ( >e/ ). 



□ 

Remark that we proved that the pair of variables (F(n re i), F(j re i)) is dynamically inde- 
pendent from the other variables of the system. The pair behaves dynamically as a pair of 
quantum oscillator variables with a frequency equal to 2 7 . 

Proposition 4 (Virial theorem). The mean square fluctuation of the relative number 
operator is proportional to the mean square fluctuation of the relative current operator, in 
particular: 

Co(F(n rel ) 2 )=(2 1 ) 2 u(FUrei) 2 ). 
Proof. This follows from the time invariance of u, i.e. u) o 8 U = 0: 

= u>(L[F(n rel )F(j rel )]) =u>(6 u [F(n rel )]F(j re i)) + u(F{n ra )~8 w [F{j rel )\) , 
and the equations of motion ( [25]) and (|26|) . □ 

Proposition 5. The mean square fluctuation of the relative number operator is given by 

u(F(n re i) 2 ) = c re / 7 coth/3 7 . 

Proof. We compute this quantity using the correlation inequalities (|T0|), rewritten in the 
form 

u(XX*) ~ n cu(XX*) ~ u{X*X) ' 1 ' 

We take for X the operator Al = FL(n re i) + i(2 7 )Fi(jVeO and then let L — ► oo, and use 
proposition 0. 

One gets 

lim u{A L A* L ) = uj(F(n rel ) 2 ) + {2 1 ) 2 u{F{j rel f) + (2 7 )c re ,, 
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and by the virial theorem (proposition [|) : 

lim oo(A L A* L ) = 2u(F(n rel ) 2 ) + (2 7 )c rei . 



L—fOO 



Analogously 



lim u{A* L A L ) = 2u(F(n rel ) 2 ) - (2 7 )c reZ . 

L — »oo 



On the other hand, 

L{A L ) = -i(2-y) 2 F L (j rel ) - (2 1 )F L (n rel ) = -(2j)A L , 



and hence 



L^oo 

lim 

L^oo 

After substitution in fl2T|) one gets 



lim w(A* L 6 u (A L )) = -4~/u(F(n rel ) 2 ) + (2 7 ) 2 c rei 
lim u)(A L 8 u {AD) = ^Co{F{n rel f) + (2 7 ) 2 c re/ . 



in 2u(F(n rel ) 2 ) -(2 7 )c re , = _ 
U 2u(F(n rel ) 2 ) + (2 7 )c re z~~ 



or alternatively 



(F(n reZ ) 2 ) = c re/7 coth/3 7 . 



□ 



This finishes the complete study of the static and dynamic properties of the canonical pair 
(F(n re i), F(j re i)) of the relative density and current fluctuations. We proved rigorously that 
for all temperatures below the condensation temperature and with non-zero condensate, 
this pair behaves like a pair of quantum harmonic oscillator variables, describing oscillations 
of the fluid from type 1 into type 2 and vice versa, yielding the typical interference pattern. 
The plasmon frequency is given by 2 7 . All this is physically clear. 

Our next and final problem is to find out what the position of the phase is in all this. We 
turn our attention now to look for a relation between the relative current fluctuation op- 
erator F(j re i) and the relative phase fluctuation operator, which we define in an analogous 
form as the total phase fluctuation operator, as follows: 



Denote its central limit by F((j) re i). 
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First, observe that one can distinguish two terms in the relative current fluctuation (|24|), 
namely the k = part and the rest: 



F L (jrel) = W-T7m(K,oKo ~ &l,oM + TTTTm Yl h l,k h ~,k ~ &!,*&+>*■ 



Denote the first term by 



F ^Urel) = 2 7 yl/2 (K,0 b -,0 ~ &-,0 & +,o), 



and its central limit by F(j® el ). Also denote by u> 9 the ground state, obtained as the 
zero-temperature limit of the equilibrium state u: 

= lim 

and cD ff the corresponding ground state for the limiting fluctuation operators observables: 

= lim 

p— >oo 

Proposition 6. We have the following relationships between the limiting fluctuation op- 
erators: 

(i) V/3 > ; p — oo included, 



F{frel) = ^(<M 

7 



^ for P = oo, 



F(jrel)=F(f rel ) = ^F((f> rel ). 

7 



Proof. As shown in |||9|], two fluctuation operators F(A),F(S) are equal in the algebra 
of fluctuation operators whenever 



(F{A-B) 2 )=0, (28) 



i.e. whenever the variance of the difference A — B of the operators vanishes. This is 
expressing in a mathematical rigorous setting, the phenomenon of coarse graining on the 
level of fluctuations. 

Therefore we calculate 



^HnfredFiM) - ^0)(F(^)F(j r J). 
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One finds, using the explicit knowledge of the state u (section |3[): 

u{F(f rel f)=u{F(<p rel f) = |^coth/? 7 



u{F(j rel )F(<p rel ))=u(F(<p rel )F(f rel )) = ^coth/? 7 , 
hence leading to the following equality, as operators: 



7 

From proposition |] and [5], it follows that 



u{F( Jrei y) = C -^coth(3 1 , 



and from proposition [|, 



lim C re l = —. 



Therefore 



This implies necessarily 



uJ 



I rel 



'([F(j rel )-F(f r 
and again the equality of the operators 

F(jrel) = F(f rel ) 

in the ground state, as a result of coarse graining. □ 

The physical interpretation of this proposition is the following. For non-zero temperatures, 
the relative current consists of two terms. One of them is j® el , which has a non-trivial 
contribution to the fluctuation of the relative current only if po > 0, i.e. whenever the 
gauge symmetry is spontaneously broken. The other term contains no more reference to 
the zero mode, in other words to the condensate. Therefore it is clear that the fluctuation 
operator F(j® el ) contains all the information of the fluctuations of what one could call 
the condensate current, or the current between the condensates interacting through the 
Josephson junction. The important equality 



F{f rel ) = #F 



rel ) 

7 
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is nothing but a rigorous translation, on the level of the fluctuations, of the popular state- 
ment: "the (superfluid, condensate) current is the gradient of the phase". The second 
statement of the proposition shows that the quantum effects on the level of the fluctua- 
tions, originating from the spontaneous symmetry breaking, are only present in the ground 
state. This of course is popular wisdom, already experienced in many models fT0|, |T4|. f22[ . 
but expressed here in a mathematically rigorous fashion for our model. 

Finally, it may come as a surprise that our results show no dependence on ip, the expectation 
value for the phase difference between the condensates (see equation (|15D). This however is 
a simple consequence of the description of the system in its mathematically simplest form, 
using the operators b± k (H), which yield a (^-independent description of the system. Indeed, 
from a mathematical point of view, the system can not be expected to behave different for 
different </?, since e.g. the eigenvalues of the Hamiltonian E k are ip- independent. 

If one is interested in the physics following from a non-zero ip, i.e. if one wants to derive 
typical Josephson currents proportional to sin ip, one needs to work with the bare operators 
ajL 2 n k . In particular, consider the relative current fluctuation operator defined by 

in stead of (|22|). 

It can easily be calculated that this operator satisfies 

lim [F L ( Jrel ),F L (j? el )]=0 

lim [F L (n rel ),F L (ff )] =ic rd smip 

5 u {F L {f rel ))=ismipF L {n rel ). 

Together with the results above, this establishes that the limiting fluctuation operator 
F Urd) is S iven b y 

F(j? el ) = F(j rel )sm<p, (29) 

where (P§|) is to be understood in terms of the equivalence between fluctuation operators 
(p8|), i.e. (p9|) follows from 

lim J(F L {j^) ~ sm<pF L ( 3rel )) 2 ) = 0. 

The variances of F(jf el ) and its dynamics, computed from proposition [5] show the explicit in- 
dependence, and its typical 'collapse and revival' properties, found in the experiments |0 . 
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